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Figure 1. Flow diagram for the Century model. 
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Figure 2. Plant production submodel flow diagram for grassland systems. 
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fine branches, large wood, and coarse roots. Wood decomposition is 
simulated in the forest model for three dead wood boxes (large wood, 
fine branches, and coarse roots). A detailed description of Century is 
presented in two papers (Parton et al., 1987, 1988) and in a user’s 
manual for the PC version. 

During the last few years, a number of simple ecosystem models 
have been developed that are amenable to regional studies over a 
long time period. Linkages (Pastor and Post, 1986) is a forest growth 
model that has been used extensively to simulate growth of different 
forest species and feedbacks between litter quality and plant growth. 
A comparison of the two models shows that soil texture, plant nitro- 
gen and lignin contents, and climatic factors are the major controls 
over nutrient cycling in both models. A major difference is that Link- 
ages uses a litter cohort approach, where the decomposition of each 
year’s litter cohort is represented separately, while in Century new 
litter is aggregated into either the metabolic or the structural litter 
component (see Figure 3). Another major difference is that Linkages 
has one SOM pool while Century has three pools. We (J. Pastor and 
W. Parton) plain on making a formal comparison between Linkages 
and Century and anticipate that the utilization of the litter cohort 
approach will be more important for forest systems since litter qual- 
ity (i.e., lignin and nitrogen content) is much more variable for forest 
species and litter types (e.g., large wood, fine wood, fine roots, and 
leaves) than for those observed in grassland ecosystems. 

The model structure used in Century is similar to that used in 
some new simple ecosystem models such as the Vegie model (Aber 
et al., 1991) and the General Ecosystem Model (GEM) (Rastetter et 
al., in press). There are a substantial number of similarities between 
the Century, Vegie, and GEM models; however, decomposition of 
plant residue and plant production are calculated using different 
approaches. The importance of model differences is not clear, and 
we plan on making a formal comparison of these models. 


Model Simplification and Parameter Determination 

Century was developed by using a simple representation of bio- 
logical processes to represent the dynamics of plant-soil systems. 
Two major model simplification processes were used: the hierarchi- 
cal approach and the conceptual procedure. 

Hierarchical Model Simplification Approach 

The hierarchical approach is based on ideas presented by Kittel 
and Coughenour (1988) and Allen and Starr (1982). This approach 
takes advantage of high-resolution models with fine time and space 
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Figure 3. Comparison of Century model structure with cohort models (from 
Schimeletal, 1991). 


scales to develop simplified relationships for coarser-scaled models. 
Model output from a fine-scaled model is transformed into an equa- 
tion that is appropriate for a coarse model. 

The hierarchical approach was used to derive the equation for the 
effect of monthly precipitation on the decomposition rate. A daily 
decomposition model and daily soil water and temperature models 
(Parton, 1978, 1984) were used to simulate daily decomposition rates 
for a 40-year time series using observed weather data to drive the 
models. These simulated daily decomposition rates were then aggre- 
gated on a monthly basis and used to develop a simplified equation 
for the effect of monthly precipitation on decomposition. A nonlinear 
least squares data-fitting procedure (Powell, 1965) was used to esti- 
mate the coefficients in the equation (Figure 4). This relationship is 
being refined in a new version of Century to include the effect of soil 
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texture on modifying the relationship shown in Figure 4. This hierar- 
chical approach provides a methodology to determine the effect of 
precipitation on decomposition for a time step (monthly) that was not 
easily derived from existing decomposition data. 

Conceptual Model Simplification Procedure 

A conceptual procedure was used to develop the overall structure 
of Century (Figure 1). This process incorporates the essential con- 
cepts of more detailed nutrient cycling models that are needed to 
simulate the dynamics of the soil system for a monthly time step. 
Simplication of the model structure is based on extensive experience 
from field and laboratory studies in establishing these concepts and 
provides the basis to extend the concepts over longer time domains 
(decades to centuries) and over greater and more diverse geographic 
regions. 

The structure of Century represents a simplification of concepts 
found in the process-oriented nutrient cycling models developed by 



TPI/PET 


Figure 4. The effect of the ratio of precipitation (ppt) to potential evapo- 
ration (pet) on monthly decomposition rates (from Parton et al., 1987). 
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Hunt (1977) and McGill et al. (1981). McGill et al.’s model (Phoenix) 
is a detailed process-oriented nitrogen cycling and soil organic mat- 
ter model (Figure 5) with a 0.2-day time step. The major structural 
simplification of the Phoenix model was accomplished by combining 
bacteria, fungi, and microbial products into a single state variable 
(active SOM). The specific effects of bacteria and fungi on the system 
were incorporated by using different microbial growth efficiencies for 
plant material decomposed at the soil surface and that decomposed 
in the soil. Surface litter is primarily decomposed by fungi, while soil 
litter is primarily decomposed by bacteria. This allowed us to sim- 
plify the model and still include the major role of bacteria and fungi. 
The conceptual simplification and aggregation of state variables in 
Century have resulted in the ability to test the potential role of 
microbial populations on SOM decomposition. In this case, the 


live roots 

V y' 

live shoots , 

\ / 


standing dead 


standing dead 

(metabolic) 

\ 

(structural) 

k . 

1 

“litter” 

J > 

“litter” 

(metabolic) 


(structural) 

sssssvk 

/ 

— , 



bacteria and 

humads 

me 

actinomycetes 

■■0M 



resistant 
soil organic 


fungi 

mauer 



1 



Figure 5. Carbon flow diagram for Phoenix soil organic matter model (from 
McGill etai, 1981). 
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model is testing a hypothesis that is based on literature data that is 
not definitive, and it could potentially be used to formulate a new 
hypothesis about the role of fungi and bacteria during decomposi- 
tion. The Century model is also now being used to test the potential 
impact of earthworms on SOM dynamics. Some of the model para- 
meters that are being modified include the decay rate for slow SOM 
and microbial growth efficiency for slow SOM. 


Parameters Used 

Most of the parameters in Century were determined by fitting the 
model to 15 long-term soil incubations where different types of plant 
material were added to the soil. A nonlinear data-frtting procedure 
(Powell, 1965) was used to determine the specific coefficients used 
in the model. A major contribution of the Century model was to 
Include the effect of the soil texture on the stabilization of SOM. A 
comparison of the observed (Sorenson, 1981) and simulated effect of 
soil texture on soil C stabilization is shown in Figure 6. It is impor- 
tant to note that model coefficients determined during the original 
model fitting procedure have not been changed. 

Model Testing and Development 

The initial version of Century (Parton et al„ 1987) was tested by 
comparing simulated long-term average plant production and steady 
state soil C and N levels with observed data from the U.S. Great 
Plains. The observed vs. simulated comparisons for soil C (Figure 
7a) and plant production (Figure 7b) show that the model did an 



Figure 6. A comparison of observed and simulated effect of soil texture 
on soil C stabilization (from Parton et al, 1987). 
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excellent job of simulating long-term grassland plant production for 
the Great Plains and did an adequate job of simulating soil C levels 
for different soil textures and abiotic environments in the region. 
One of the major parts of the model testing procedure was an analy- 
sis of observed regional plant production data (Sala et al., 1988) and 
soil C and N data (Burke et al., 1989) to develop a validation data 
set, since an adequate regional data set was not otherwise available 
for testing the model. 

This regional soil data base for the Great Plains (Burke et al., 
1989) was used to test and validate the ability of Century to simulate 
regional patterns of soil C and N for grassland soils (Parton et al., 
1987). This data base is now being used to test new concepts about 
the formation of passive soil C (old soil C) and to test the ability of 
Century to simulate C and N losses due to cultivation. Recent data 
indicate that the slow SOM pool is not the primary source of passive 
soil C (the assumption made during the initial model development) 
and that the amount of passive soil C is higher for clay soils. The 
inability of Century to predict SOM dynamics for forest soils and to 
predict observed differences in C-to-N ratios of SOM for different soil 
textures (i.e., 13 to 14 for silty soils and 10 to 1 1 for clayey soils) has 
also helped to formulate these concepts. In the recent formulation, 
most of the passive soil C is derived from the turnover of active SOM, 
and passive SOM increases with clay content. By incorporating these 
concepts, the model correctly predicts observed differences in soil C- 
to-N ratios and predicts that the fraction of soil C in the passive pool 
will range from 30% for silty soils to 50% for clay soils. The regional 
soil data base is being used to test how well the model simulates soil 
textural differences in C-to-N ratios and steady state soil C and N 
levels. This demonstrates the importance of using regional data 
bases to test model performance and to help reformulate the model. 
Testing of general ecosystem models with a large number of sites is 
an essential part of the model development and testing procedure. 

The original version of the model was tested by comparing simu- 
lated plant production and soil C and N levels with extensive data 
on plant production and soil organic matter (C and N) from the U.S. 
Great Plains (Parton et al., 1987). We are presently testing the model 
by using grassland plant production and soil organic matter data 
from 15 sites around the world. The testing of the model using data 
from many different grassland ecosystems has been an important 
part of the model development process and has suggested model 
changes which were necessaiy to formulate a truly general grass- 
land ecosystem model. 

During the last two years we have developed a more mechanistic 
version of the grassland plant production submodel with the objec- 
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tive of being able to simulate year-to-year variations in annual plant 
production and to simulate the seasonal dynamics of live plant bio- 
mass. The new model calculates maximum monthly plant produc- 
tion as a function of soil temperature, monthly precipitation, and 
the ratio of live shoots to dead shoots. The model also includes the 
effects of grazing and fire on plant root-to-shoot ratio and N content 
of the plant material (Holland et al., in press; Ojima et al., 1990). 
This new model was developed using long-term (1928-89) plant pro- 
duction data from a tallgrass prairie in eastern Kansas (Towne and 
Owensby, 1984; Owensby, unpublished data) and a seven-year time 
series of plant production data from a shortgrass prairie site in 
Colorado (Dodd and Lauenroth, 1978). Researchers measured plant 
production from 1970 to 1976 for control, irrigated (only), fertilized 
(only), and irrigated plus fertilized sites. Simulated annual plant 
production for the different treatments in the shortgrass prairie (Fig- 
ure 8a) and a comparison of observed and simulated production 
(Figure 8b) show that the model reasonably simulates plant produc- 
tion for the different treatments (r 2 = 0.85). A comparison of 
observed and simulated annual aboveground plant production data 
from 1930 to 1968 for the tallgrass prairie site (Figure 9a) shows 
that the model has an observed and simulated r 2 of 0.60. Note that 
the observed plant production for these years is underestimated 
since the biomass measurements were made by clipping the bio- 
mass at a height of 5 cm, thereby omitting a portion of the grass 
biomass. Figure 9b shows the comparison of observed and simu- 
lated biomass data from 1970 to 1989 for the late-sprlng-bumed 
and unbumed sites. The model correctly estimated plant production 
for the very dry years (1981 and 1989) and correctly predicted the 
mean difference between the burned and unbumed sites. Year-to- 
year variability was not as well predicted by the model; however, it is 
important to note that the standard deviation associated with the 
production estimates ranges from 100 to 200 g/m 2 /yr and makes it 
difficult to assess the significance of differences between modeled 
and observed data. In summary. Century reasonably simulated the 
Impacts of fertilization, irrigation, and burning on plant production 
over an extended period of time, but further refinements may be 
needed to simulate smaller year-to-year variations in production. 

The ability of the model to simulate seasonal biomass dynamics 
for live and standing dead biomass Is now being tested by compar- 
ing the model results to monthly plant biomass data from four tropi- 
cal grassland sites (Mexico, Kenya, Ivory Coast, and Thailand) and 
temperate grassland sites in the former USSR and Ireland. The 
results for the tropical sites are very promising and show that the 
model can correctly simulate the seasonal dynamics of live and dead 
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biomass and respond appropriately to fire events. The turnover rate 
of live roots, the rate of transfer from standing dead biomass to lit- 
ter, and the root-to-shoot ratio are parameters that vary consider- 
ably among the sites. These differences are caused by plant-specific 
differences and abiotic factors such as snow. The observed and sim- 
ulated live biomass for the Lamto site in Ivory Coast agree reason- 
ably well (Figure 10). The results for the Russian grassland sites 
indicate that the model underestimated plant production and soil N 
inputs to these systems and suggest that the equations for the effect 
of rainfall on plant production and soil N inputs need modification. 
Such testing of the model at sites with diverse climatic conditions is 
an important part of model testing and validation and has led to the 
development of a more accurate and general grassland model. 

Regional Modeling 

The Century model has been used to simulate regional patterns of 
grassland biogeochemistry (Burke et al., 1990; Parton et al., 1989). 
A hierarchical approach is used to simulate the regional ecosystem 
dynamics (Figure 11). Century is a patch model which represents 
plant production as an aggregate variable from a mixture of plant 
species. The external driving variables for the model are soil texture 
and monthly climatic data. In scaling up from patch to topose- 
quence in a landscape or to a physiographic unit, the diversity of 
soil types tends to increase. In making a regional simulation with 
Century, soil variability is accounted for by identifying key soil types 



Figure 1 0. Comparison of observed and simulated live biomass for a 
burned grassland site in West Africa . 
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and parameterizing Century for each individual type. The results 
can then be aggregated as an areal average weighted mean. In addi- 
tion, the range in climatic factors increases with the size of a region. 

Overlaying regional soil variability on climatic variability within a 
region defines unique sets of driving variables for Century. This is 
facilitated by a geographic information system (GIS) that contains 
information on the spatial patterns of soil texture and monthly cli- 
matic data. A regional grassland simulation for northeastern Colo- 
rado (Burke et al., 1990) required 160 model runs to represent the 
unique combinations of soil texture and climate. 

Across larger areas, dominant plant species changes may result 
in modifications to certain key ecosystem characteristics, such as 
the root-to-shoot ratio and plant lignin and nutrient content. For 
the Great Plains region we have developed equations which predict 
how the root-to-shoot ratio and plant lignin and nutrient contents 
change as a function of the annual precipitation. These equations 
work because plant communities change as a function of annual 
precipitation. 

Simulated annual patterns for aboveground production, soil 
organic C, net N mineralization, and N gas flux (N 2 0, NO x , and NH 3 ) 
in northeastern Colorado (Figure 12) show that soil organic C is pri- 
marily controlled by soil texture, while regional patterns of the other 
variables correspond to annual precipitation (Burke et al., 1990). The 
effect of spatial resolution of input data indicates that regional esti- 
mates of the average plant production and N gas flux are relatively 
insensitive to changes in spatial resolution of inputs from 2.5 to 
4000 km 2 . In contrast, average soil C levels are substantially under- 
estimated (14%) when the spatial resolution of inputs is decreased 
from 1500 to 4000 km 2 . These results suggest that aggregation 
errors are noticeable for soil C levels when the spatial grid is coarse 
because of the nonlinear impact of soil texture on soil C levels. 

Regional patterns for grassland ecosystem properties have been 
simulated for the North American Central Grasslands (Parton et al., 
1989; Schimel et al., in press) based on regional variation in cli- 
mate. While results are limited by the lack of detailed soil and land 
use data, the regional pattern of net primary production compares 
well with that of the annual integral of satellite-based vegetation 
index data (Schimel et al., 1991). At present, development of a 
regional GIS data base that includes climatic, land use, and soil 
data for the North American Central Grasslands is ongoing, and the 
data base will be used to more accurately simulate patterns of soil C 
and plant production. 

Sensitivity of ecosystem processes and of fluxes of C and N to 
changes in climate and atmospheric C0 2 level has been evaluated 
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Figure 12. Simulated output from the Century ecosystem model for (a) net 
primary production in g C/mVyr, (b) soil organic carbon in g C/m to 20 cm 
depth, (c) net annual 1 V mineralization in g N/m^/yr to 20 cm depth, and 
(d) net annual N gas flux in northeastern Colorado in g N/m /yr. Lines 
overlaid represent annual precipitation contours (a, c, d) and soil texture 
classes (b) (from Burke et al, 1990). 


for Century. Results from simulating the effects of a GCM C0 2 -dou- 
bling climate scenario (Hansen et al., 1988) suggest that plant pro- 
duction will generally increase in the northern portion of the North 
American Central Grasslands and decrease in the southern portion 
largely in response to precipitation (Schimel et al., 1990; Kittel et 
al in preparation). In addition, throughout the region there will be 
a net release of C from the soil and an increase in fluxes of nitroge- 
nous gases (N 2 . N 2 0. and NH 3 ). 


Conclusions 

There has been substantial progress in the development of simpli- 
fied ecosystem models which have the potential to be incorporated 
into earth system models and to address regional issues regarding 
changes in land use or climatic conditions. An essential part of the 





‘liWiam J. ‘Parton ef at. 


299 


development of these models Is detailed comparison of models with 
data from a large number of sites having different climate and soil 
characteristics. These comparisons are critical for testing this class 
of models and have greatly enhanced our understanding of the com- 
plex interactions, leading to substantial improvements in their struc- 
ture. These models are capable of simulating the impact of climatic 
change on different ecosystems and can provide insight to potential 
changes in carbon (e.g., C0 2 ) and nutrient (e.g., N 2 0) fluxes from 
ecosystems that, as radiatively active gases, play a role in global cli- 
mate. A significant problem with including simple ecosystem models 
in earth system models is that the spatial scale of the atmospheric 
GCM models is still too coarse to incorporate the spatial heterogene- 
ity in surface climate and land surface properties (e.g., topography) 
needed to represent natural ecosystems. This is illustrated by the 
fact that the southern part of the Great Plains, which includes short- 
grass, mid-grass, and tallgrass prairie ecosystems; montane systems; 
and agroecosystems that range from wheat-fallow to continuous corn 
systems, is covered by a single grid with 7.8° latitude x 10° longitude 
in the GCM at the Goddard Institute for Space Studies. Even with a 
very fine GCM grid network, it is anticipated that a need to consider 
subgrid diversity in agricultural land use and natural ecosystems 
must be explicitly addressed. This level of resolution will be neces- 
sary to represent the potential interactions between the atmosphere 
and the biosphere (Pielke and Avissar, in press; Avissar and Pielke, 
1989; Avissar and Verstraete, 1990). 
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